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FOREWORD 


In  this  report,  a  procedure  is  outlined  which  allows  the  generation  of  an 
S-N  curve  with  some  statistical  value  without  resorting  to  an  extremely 
large  data  base. 


This  report  was  prepared  in  the  Mechanics  and  Surface  Inte..  ions 
Branch  (AFML/MBM) ,  Nonmetallic  Division,  Air  Force  Materials  Laboratory, 
Wright-Patterson  AFB,  Ohio.  The  work  was  performed  under  Project  2419, 
"Nonmetallic  Structural  Materials,"  Task  No.  241903,  "Composite  Materials 
and  Mechanics  Technology."  The  time  period  covered  by  this  effort  was 
from  1  October  1978  to  15  April  1979.  James  M.  Whitney  (AFML/MBM)  was  the 
laboratory  project  engineer. 

The  author  wishes  to  acknowledge  Professor  W.  Park  of  Wright  State 
University  for  helpful  suggestions  concerning  statistical  analysis  and 
G.  Sendeckyj  of  the  Air  Force  Flight  Dynamics  Laboratory  for  useful  dis¬ 
cussions  concerning  "strength  degradation"  and  "wearout"  models. 
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SECTION  I 
INTRODUCTION 

The  classical  S-N  curve  has  been  the  primary  method  of  characterizing 
the  fatigue  behavior  of  fiber  reinforced  composites.  This  method  usually 
consists  of  determining  the  number  of  cycles  to  failure  for  a  number  of 
maximum  stress  ranges  associated  with  a  particular  load  history  (often 
constant  amplitude  tension-tension  loading) .  The  resulting  S-N  curve 
yields  an  estimate  of  the  mean  time-to-failure  as  a  function  of  maximum 
stress  range.  Such  a  procedure,  however,  fails  to  account  for  the  large 
variation  in  the  time-to-failure  at  a  given  maximum  stress  level.  Fatigue 
data  with  statistical  significance  requires  a  large  number  of  replicates 
at  a  given  maximum  stress  level  in  order  to  measure  the  distribution  of 
t ime-to-f ailure . 

In  this  report  a  procedure  which  allows  the  generation  of  an  S-N 
curve  with  some  statistical  value  without  resorting  to  an  extremely 
large  data  base  is  explored  in  detail.  The  "wearout"  or  "strength 
degradation"  model  approach  (References  1,  2,  3)  provides  one  means  for 
accomplishing  this.  Such  an  approach,  however,  involves  the  assumption  of 
a  direct  relationship  between  static  strength  distribution,  residual  strength 
distribution  after  time  under  a  specified  load  history,  and  distribution  of 
time-to-failure  at  a  maximum  stress  level.  These  types  of  models  have  to  be 
carefully  defined  in  terms  of  the  load  history  to  be  applied  to  the 
material.  For  example,  if  the  load  history  is  tension/compression  then 
the  concept  of  strength  degradation  must  consider  both  residual  tension 
and  residual  compression,  along  with  possible  competing  failure  modes. 

The  alternative  approach  to  S-N  curve  characterization  does  not  require 
any  assumptions  relating  fatigue  life  to  residual  strength.  As  shown  in 
the  Appendix,  the  procedure  does  measure  parameters  which  are  completely 
compatible  with  the  "wearout"  model  approach  without  any  residual  strength 


measurements. 
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SECTION  II 

S-N  CURVE  CHARACTERIZATION 


An  alternative  to  the  "wearout"  model  approach  for  fatigue  charac¬ 
terization  of  composite  materials  has  been  proposed  by  Hahn 'and  Kim 
(Reference  4).  Their  approach  involves  two  basic  assumptions:  (1)  a 
classical  power  law  representation  of  the  S-N  curve  (Reference  5),  and 
(2)  a  two  parameter  Weibull  distribution  of  time-to-failure.  In  mathe¬ 
matical  form  these  assumptions  become 

CNSb  =  1  (1) 

where  S  is  the  stress  range,  N  is  the  number  of  cycles  to  failure,  and 
b  and  C  are  material  constants.  In  addition. 


R  (N)  =  exp 


(2) 


where  R  (N)  denotes  the  reliability  of  N  (probability  of  survival) ,  Nq 
is  the  characteristic  time-to-failure  (location  parameter),  and  is 
the  fatigue  shape  parameter.  As  illustrated  in  the  Appendix,  Equations 
1  and  2  can  also  be  derived  from  the  "wearout"  model  approach  with  . 
being  related  to  the  shape  parameter  for  static  strength  (Reference  3) . 
Equation  1  can  be  written  in  the  form 

- 1  -b 

N  =  C  S  (3) 

substituting  Equation  3  into  Equation  2  and  solving  for  S  yields 


S 


«f  *  f(S) 


where 


When  N  =  N  ,  -In  R  (N  ) 
o  o 


1  ’ 

K  =  C  b 
and  Equation  4 


reduces  to 


S(N  )=KN 
'  o'  o 


(4) 

(5) 

(6) 


2 


AFML-TR-79-4111 


A  plot  of  log  S  versus  log  N^  produces  a  straight  line  with  slope  **  1/b 
and  a  y  intercept  of  log  K.  Thus,  a  measure  of  the  distribution  of  timer* 
to-failure  at  various  stress  ranges  in  conjunction  with  Equation  2  allows 
ci^  to  be  determined  along  with  a  set  of  values  of  Nq,  each  corresponding 

to  a  value  of  S.  Equation  4  can  then  be  utilized  to  produce  an  S-N  curve 

••  o 

for  any  desired  reliability,  R(N) . 

However,  from  a  practical  standpoint,  one  is  more  interested  in 
obtaining  an  S-N  curve  for  any  desired  level  of  reliability  rather  than 
an  S-N^  curve.  Writing  Equation  6  in  the  form 

N0=C_1s“b  (7) 

and  substituting  into  Equation  2  yields 


R(N)  -  exp 


a, 


(8) 


Solving  for  S  leads  to  the  following  S-N  relationship  for  any  desired  level 
of  reliability 


S  =K 


(9) 


3 


il 


>  .' 


iy*., 


.  •  _•; 


* 


AFML-TR-79-4111 


SECTION  III 


DATA  REDUCTION  PROCEDURE 


The  data  reduction  procedure  consists  of : 


(a)  Fitting  the  time-to-failure  data  at  each  stress  range  to  a  two 
parameter  Weibull  distribution. 


(b)  Use  a  data  pooling  scheme  to  determine  the  fatigue  shape 


parameter  a^. 


(c)  Fit  log  S  versus  log  N^  data  to  a  straight  line  for  the  deter¬ 


mination  of  b  and  K. 


Let  ra  be  the  number  of  stress  ranges  tested  and  n  the  number  of 


specimens  tested  at  the  i-th  stress  range,  S^,  which  leads  to  the  data 


N.  (N  ,  N.  ,  • « • »  N.  ),  l-l,  2,  •  •  • »  m 

ill  lb  I 


Each  stress  range  is  fit  to  the  two  parameter  Weibull  distribution 


R(N. )  =  exp 


[-&P 


A  number  of  procedures  can  be  utilized  for  determining  and  Nq^. 


One  of  the  methods  preferred  by  statisticians  is  the  maximum  likelihood 
estimator  (MLE)  which  is  of  the  form  (Reference  6) 


EN..  «fi  In  N.. 
l3  l3 


MX 

M 


« 

H 


l» 


v 
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(13) 


A  *  - 

where  ct^  and  NQi  denote  estimated  values  of  and  NQi,  respectively. 
Equation  12  has  only  one  real  positive  root.  As  a  result,  an  iterative 
scheme  can  be  utilized  until  a  value  of  ot  is  obtained  to  any  desired 
number  of  decimal  places.  The  resulting  value  of  obtained  from  the 
iterative  scheme  can  then  be  used  in  conjunction  with  Equation  13  to 

A 

obtain  N  , . 

oi 


Since  fatigue  shape  parameters  are  estimated  based  on  a  particular 
sample  size,  it  is  anticipated  that  each  value  of  would  produce  a 
different  value  of  even  though  may  be  independent  of  stress  range. 

A  two  sample  test  (Reference  7)  is  available,  however,  which  allows  for 
testing  the  equality  of  shape  parameters  in  two  parameter  Weibull  distri¬ 
butions  with  unknown  scale  parameters.  The  approach  is  based  on  MLE  and 
the  results  depend  on  sample  size  and  confidence  level  desired.  Let 

a  and  a  .  be  the  maximum  and  minimum  values  obtained  for  oL , .  For 
max  min  f  i 

the  information  tabulated  in  Reference  7,  it  is  required  that  a  and 

’  max 

a  .  be  associated  with  equal  sample  sizes,  n.  If  a,  and  a.  .  are 
min  n  r  »  f  max  f  min 

from  the  same  distribution,  then  it  is  expected  that  (Reference  7) 


A 

0(y  ,  n),  0>1  (14) 

af  min 

for  a  given  confidence  level,  a,  and  sample  size,  n.  Values  of  3  are 
shown  in  Table  1  for  various  sample  sizes  corresponding  to  a  confidence 
level  of  0.98.  This  data  is  taken  from  Reference  7.  The  large  values 
of  A  associated  with  small  sample  sizes  suggest  that  significant  variations 
in  a  ^  are  likely  to  be  encountered  with  small  data  sets  taken  from  the 
same  population. 
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Let  us  now  assume  that  is  independent  of  stress  range.  Then  a 
data  pooling  technique  must  be  utilized  in  order  to  determine  a  single 
value  of  for  all  S^.  Various  approaches  for  obtaining  a  pooled  value 
of  can  be  found  in  the  literature.  The  approach  used  in  the  present 
work  has  been  investigated  by  Lemon  (Reference  8).  This  procedure  utilizes 
the  normalized  data  set 


N 


(15) 

(16) 


Thus,  each  set  of  data  at  a  given  stress  range  is  normalized  by  the 
estimated  characteristic  time-to-failure  and  the  results  fit  to  the  pooled 
two  parameter  Weibull  di stribution 


This  procedure  has  the  advantage  of  obtaining  a  large  data  base  for 
determining  by  using  a  few  replicates  for  a  number  of  values  of  S^. 
In  general,  for  equal  accuracy  fewer  specimens  are  needed  to  determine 
tne  location  parameter  than  shape  parameter 


For  the  pooled  Weibull  distribution,  Equation  17,  the  MLS  relation¬ 
ships  tane  the  form 


6 


AFML-TR-  7  9-4  111 


where  and  Xq  are  estimated  values  of  and  Xq,  respectively,  and 


M  = 


■Ei 


For  a  perfect  fit  to  the  data  pooling  scheme,  the  location  parameter, 

A 

XQ,  should  be  unity.  The  value  of  can  be  adjusted  to  produce  an 
exact  value  or  unity  for  Xq.  In  particular, 

N  .  =  X  N  .  (21) 

Ol  O  Ol  '  ' 

where  N  .  denotes  estimated  values  of  N  associated  with  the  adjusted 
oi  o 

two  parameter  Weibull  distribution 


R  (X)  =  exp  (-X°S:) 


The  slope  of  the  S-N  curve,  1/b,  and  the  y-intercept,  K,  can  be 
determined  by  fitting  log  versus  log  NQi  to  a  straignt  line.  With  K, 
b,  and  now  determined.  Equation  9  can  be  used  to  produce  an  S-N  curve 
of  any  desired  reliability. 


It  should  be  noted  that  MLE  is  asymptotically  unbiased,  i.e.,  it  is 

a  biased  estimator  for  small  sample  sizes  (Reference  6) .  Unbiasing 

factors  are  tabulated  in  Reference  9.  These  factors  are  less  than  unity 

as  MLE  always  tends  to  overestimate  the  shape  parameter.  Confidence 

intervals  for  both  the  shape  parameter  and  location  parameter  as  a  function 

of  sample  size  have  also  been  established  (Reference  9).  If  conservative 

estimates  are  desired  for  R  (N) ,  then  a  lower  bound  value  of  can  be 

utilized.  This  value  of  can  be  used  in  conjunction  with  Equation  19 

to  determine  X  . 

o 


•'M 

.38 


7  $ 
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SECTION  IV 

CENSORING  PROCEDURES 

In  the  case  of  high  cycle  fatigue  the  time-to-failure  may  become 
unacceptably  long.  This  difficulty  can  be  overcome  by  raising  the  stress 
range  so  that  fatigue  failures  are  produced  within  a  reasonable  number  of 
cycles.  For  filament  dominated  laminates,  which  tend  to  have  a  very  flat 
S-N  curve,  stress  levels  may  have  to  be  raised  to  an  unacceptable  level 
to  produce  a  reasonable  time-to-failure  for  all  specimens  tested.  In 
particular,  it  is  undesirable  to  raise  the  fatigue  stress  level  to  such 
a  degree  that  it  significantly  overlaps  the  static  strength  distribution. 

In  such  cases  the  probability  of  a  first  cycle  failure  is  significant. 

If  censoring  techniques  are  applied  to  data  reduction  procedures, 

fatigue  failures  are  not  required  of  all  specimens.  For  the  data  reduction 

scheme  outlined  in  the  present  work.  Type  I  censoring  seems  to  be  the 

most  desirable  in  terms  of  yielding  the  most  information.  In  the  case  of 

Type  I  censoring  the  fatigue  test  is  terminated  at  a  prs-determined  time 
£ 

(e.g.  10  cycles)  even  though  all  specimens  have  not  failed.  The  MLE 
equations  for  Type  I  censoring  are  of  the  form  (Reference  6) 


a  A 

EN„  fi  In  N..  +  (n.  -  r.)  R.  ®fi  In  R. 

ij  ij  l  l  l  l 

H _ _ 

ri  * 

LN„  fi  +  (n.  -  r.)  R.  fi 

ij  i  i'  l 

j=l 


r. 

_L_  \  ’  In  N.,  - —1—  =0 

ri  ^  *fl 


(23) 
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A 

a 

N..  fi 


(V 


A 

r.)  R.afi 

i  1 


(24) 


where  n^  now  denotes  the  total  number  of  specimens  tested  at  S^,  is 
the  number  of  fatigue  failures  at  S^,  and  is  the  number  of  cycles  at 
which  the  test  is  terminated. 


The  data  pooling  procedure  is  now  analogous  to  "progressive  censoring" 
In  which  a  number  of  samples  are  removed  at  pre-determined  time  intervals 
throughout  the  duration  of  the  test.  The  MLE  associated  with  the  data 
pooling  procedure  in  conjunction  with  censored  samples  becomes 


m  r. 


m 


EExii“'k  xu  +£"vri,Yi““n',i 

i=i  j=i 


i=l 


m 


m  r. 

LEX" +  E  ,ni  •  v  Yi 

*i=l  j=l 


(25) 


i=l 


where 
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SECTION  V 
EXAMPLE  DATA 


Consider  the  tension-tension  (T-T)  fatigue  data  presented  by  Ryder 

*> 

and  Walker  (Reference  10)  on  quasi-isotropic  T300/934  graphite-epoxy 
laminates  with  the  stacking  geometry  [45/90/-45/90-45/0/45/0]  (laminate 

D 

2) .  Twenty  replicates  at  three  stress  levels  were  used  to  characterize 

the  scatter  in  time-to-failure.  The  effect  of  sample  size  on  the  Weibull 

parameters  can  be  estimated  by  using  a  table  of  random  numbers  to  select 

sample  sizes  of  5,  10.  15,  and  20  for  the  three  stress  levels  and  calculating 

the  resulting  Weibull  parameters.  The  results  are  shown  in  Table  2  along 

with  the  pooled  parameters  X  and  a..  The  Weibull  parameters  are  determined 

o  „  f 

by  KLE.  Note  that  the  values  of  are  not  radically  effected  by  sample 

size,  while  is  very  sensitive  to  n.  The  same  t  end  is  noted  for  the 

pooled  parameters  X  and  a..  It  should  also  be  noted  that  the  data  for 
o  f 

n=2G  satisfies  the  criterion  of  Equation  14  for  data  pooling. 


Comparison  between  the  normalized  data  and  the  Weibull  distribution 
obtained  from  the  pooling  procedure  shown  in  Figure  1.  Three  additional 
stress  levels  with  3  replicates  each  are  added  to  the  data  pooling  process, 
yielding  a  total  of  69  pooled  date  points  in  Figure  1.  Comparisons  between 
data  and  Weibull  distributions  obtained  from  the  pooling  procedure  are 
shown  in  Figures  2-4  for  the  three  stres-  levels  with  20  replicates.  Tn 
Figures  1-4,  data  points  are  converted  t3  probabilities  of  survival  from 
the  Median  Rank  (MR)  defined  as 


MR 


^  j  -  0.3 
n  +  0.4 


(29) 


where  j  is  the  survival  order  number  (data  listed  'n  decreasing  order  of 
time-to-failure)  end  n  is  the  total  number  of  samp..es  tested.  Equation  29 
is  an  estimate  of  the  reliability  function  R.  When  probability  of  failure 
is  desired  Equation  29  ie  applied  with  the  data  listed  In  increasing  t#me- 
to- failure. 


11 
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While  the  pooled  data  in  Figure  1  compares  favorably  with  the  estimated 
two  parameter  Weibull  distribution,  similar  data  correlation  for  the  indi¬ 
vidual  stress  levels  in  Figures  2-4  are  much  less  favorable.  This  is  to  be 
anticipated  as  small  data  samples  in  conjunction  with  large  scatter  creates 
difficulty  in  fitting  the  data  to  any  reasonable  distribution  function. 

These  results  are  also  a  good  illustration  of  the  desirability  of  data 
pooling  for  obtaining  a  larger  sample  size  for  estimating  shape  parameters. 
The  characteristic  S-N  curve  resulting  from  the  data  reduction  scheme  is 
shown  in  Figure  5  along  with  a  95  percent  survivability  curve  calculated 
from  Equation  9.  The  numbers  in  parentheses  correspond  to  the  sample  size 
associated  with  a  specific  stress  level.  Scatter  bands  on  time-to- failure 
are  also  shown.  The  solid  dots  correspond  to  fatigue  failures  outside  the 
95  percent  survivabi-i  ity  line.  Tlie  pooled  Weibull  parameters  along  with 
the  S-N  curve  parameters  are  listed  in  Table  3. 

The  S-N  characterization  procedure  has  also  been  applied  to  other 
data  available  in  the  literature  and  the  results  are  shown  in  Figures  6-8. 

As  in  Figure  5,  sample  size  at  each  stress  level,  scatter  bands,  and 
fatigue  failures  outside  of  the  93  percent  survivability  line  are  shown. 


In  Figure  6  tension-compression  (T-w)  constant  amplitude  fatigue  data 
from  laminate  2  is  illustrated  in  Reference  10.  Data  in  Figure  7  is  also 
from  the  work  of  Ryder  and  Walker  (Reference  11) .  The  same  graphite-epoxy 
material  as  in  laminate  2  is  utilized  for  this  T-T  data  with  the  stacking 
geometry  [Q/+b5lQ^I -1*5/ Q^I-\-U5lQ ^1-^5/ 0"\ This  composite  is  designated 
laminate  1  and  presents  a  very  interesting  example  for  the  recommended 
fatigue  characterization  procedure.  Because  of  the  flat  nature  of  the 
S-N  curve,  fatigue  failures  were  difficult  to  induce.  For  the  split 
numbers  in  parentheses,  the  first  number  designates  the  total  number  of 
samples  tested  at  that  stress  level  and  the  second  number  refers  to  the 
number  of  fatigue  failures  recorded.  Thus,  this  data  had  to  be  reduced 
utilizing  the  censoring  procedure  outlined.  All  tests  were  terminated 
after  10  cycles.  It  is  interesting  to  note  the  extremely  low  value  of 
associated  with  this  laminate.  This  behavior  has  been  previously  noted 
for  filament  dominated  laminates  (Reference  12)  (laminates  containing  a 
large  percentage  of  0-degree  plies  relative  to  the  load  direction) .  The 
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straight  line  fit  of  the  S-N  curve  ia  not  as  good  for  this  laminate  as 
for  laminate  2  in  Figures  5  and  6.  The  low  number  of  fatigue  failures 
along  with  the  large  scatter  are  likely  reasons  for  the  poorer  fit. 

The  data  in  Figure  8  is  taken  from  the  work  of  Yang  (Reference  13) . 

This  constant  amplitude  T-T  data  are  for  a  [+45]  T300/5208  graphite-epoxy 

6  b 

laminate.  This  orientation  was  chosen  because  it  induces  significant 
shear  stress  relative  to  the  fiber  direction  in  each  ply.  Since  there 
are  no  0-degree  plies  relative  to  the  load  direction,  this  laminate  is 
referred  to  as  matrix  dominated.  This  composite  is  designated  laminate  3. 

The  laminate  numbers  have  been  assigned  in  descending  order  of  filament 
dominance.  It  is  interesting  to  note  the  relatively  high  value  of 
associated  with  laminate  3  compared  to  laminates  1  and  3.  Because  of  the 
higher  value  of  the  shape  parameter,  the  amount  of  shift  of  the  95  percent 
survivability  line  down  from  the  characteristic  life  line  s  much  less 
compared  to  laminates  1  and  2  in  Figures  5-7. 

To  dec  astrate  the  basic  difference  in  S-N  curve  behavior  for  constant 
amplitude  T-T  loading,  the  95  percent  survivability  lines  for  the  three 
laminates  under  consideration  are  shown  in  Figure  8.  Graphite  fibers  are 
essentially  fatigue  insensitive.  Thus,  filament  dominated  laminates  tend 
to  nave  a  flatter  S-N  curve. 

Pooled  Weibull  parameters  along  with  the  S-N  curve  parameters  are 
shown  in  Table  3  for  all  laminates  and  loading  conditions  under  consideration. 
Note  that  none  of  the  values  for  have  been  corrected  for  the  bias  of 

MLE.  For  these  pooled  sample  sizes  the  bias  correction  is  very  small. 

All  stress  ranges  in  Figures  5-9  are  normalized  by  the  characteristic 
static  strength  S^  obtained  by  fitting  the  static  strength  data  to  a  two 
parameter  Weibull  distribution. 
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SECTION  VI 


CONCLUSIONS 


A  procedure  has  been  outlined  for  characterizing  the  S-N  behavior 
of  composite  laminates  with  some  degree  of  statistical  reliability.  The 
data  pooling  procedure  used  in  conjunction  with  the  data  reduction  scheme 
offers  the  advantage  of  a  large  data  base  for  determining  a  fatigue  shape 
parameter  without  requiring  large  sample  sizes  at  each  stress  range 
considered.  Use  of  censoring  techniques  provides  data  reduction  procedures 
without  requiring  all  specimens  to  produce  fatigue  failures.  Data  presented 
shows  the  proposed  characterization  scheme  to  be  promising. 

If  the  proposed  procedure  is  to  be  utilized,  the  fatigue  experiments 
should  be  planned  accordingly.  In  particular,  it  would  be  desirable  to 
use  the  same  sample  size  for  all  stress  ranges  tested.  If  60  samples  are 
desired  for  determining  a  pooled  shape  parameter,  10  replicates  at  6 
different  stress  ranges  is  preferable  to  20  replicates  at  3  different 
stress  ranges.  In  the  former  case,  6  data  points  are  available  for 
determining  the  characteristic  S-N  curve,  rather  than  3  as  would  be 
provided  by  the  latter  case. 

In  the  data  reduction  procedure,  use  of  MLE  is  strongly  recommended 
for  determining  Weibull  parameters.  Confidence  intervals  and  other 
statistical  tools  based  on  MLE  have  been  well  established.  Furthermore, 
the  data  pooling  technique  for  censored  data  is  based  on  MLE  in  conjunction 
with  the  concept  of  "progressive  censoring." 

It  should  also  be  noted  that  the  procedure  outlined  can  easily  be 
revised  to  include  S-N  relationships  other  than  the  power  law  described 
by  Equation  1. 
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APPENDIX 


For  simplicity,  consider  the  case  of  constant  amplitude  tension- 
tension  loading.  If  one  assumes  a  crack  growth  model  as  the  basi£  for 
a  "strength  degradation"  model,  the  governing  equation  for  residual 
strength  takes  the  form 

.Jbjl-c) 

dF  =  _  _AS  F  (30) 

dn  c 


where  F  is  the  strength  after  n  cycles  -?nJ  A,  b,  and  c  are  constants. 
Thus,  there  are  three  parameters  which  characterize  residual  strength 
degradation  in  this  model.  For  the  assumption  that  the  flow  growth  is 
driven  by  a  square  loot  singularity  at  the  crack  tip  (Reference.  1) 

b  =  2  (1-c)  (31) 

and  the  number  of  parameters  is  reduced  to  two. 


Separating  variables  in  Equation  30  yields 


F  o 

o 


(32) 


Performing  the  integration  leads  to  the  residual  strength  equation 

F  C  =  FC(n)  +  ASbn  (33) 

o 

where  Fq  is  the  static  strength.  Let  us  assume  that  the  static  strength 
is  described  by  the  two  parameter  Weibull  distribution 


R(Fq)  =  exp 


(34) 


where  Sq  is  the  location  parameter  of  the  static  strength  distribution 

and  a  is  the  static  strength  shape  parameter.  Substituting  Equation  33 
o 

into  Equation  34  yields 


R(F)  =  exp 


(35) 
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Denoting  the  number  of  cycles  to  fatigue  failure  by  N^,  it  is  assumed  that 
failure  actually  occurs  when  the  stress  reaches  its  maximum  value  during 
the  fatigue  cycle.  In  terms  of  residual  strength  this  simply  states  that 
the  residual  strength  is  reduced  to  S  one  cycle  prior  to  failure,  as 
failure  occurs  on  the  next  cycle  when  loaded  to  S.  Stated  mathematically, 

r-,  _  o  „  =  w  =  N  (36) 


F  (n)  =  S  ,  n  =  Nf-1  =N 


and  Equation  33  becomes 


c  c  b  (37) 

F  =  S  +  AS  N 
o 

The  case  N  =  0  corresponds  to  a  one  cycle  failure  and  Equation  37  becomes 

S  =  Fo 

which  simply  implies  that  a  first  cycle  failure  is  a  special  case  where 
the  maximum  stress  in  the  fatigue  cycle  corresponds  to  the  static  strength. 

Using  Equation  37  in  conjunction  with  Equation  34  leads  to  the  following 
reliability  function  for  time-to-failure 


R(N)  =  exp 


- 1  c-b 
N  +  A  S 

.  -  J.  -b  c 
A  S  S  , 

L.  O  j 


where  is  the  fatigue  shape  parameter 

°f  =0b/c 


Equation  39  can  be  written  in  the  form  of  a  three  parameter  Weibull 


distribution 


R  (N)  =  exp 


rN  -  L  \ 


N  -  L 
o  o 


n 


where 


N  =S  CA' 
o  o 


vT(t)1 


t  » - 1  _c-b 
L  =  -A  S 
o 


m* 


V 

I 


Ws» 
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The  three  parameter  Welbull  distribution  is  simply  a  shifted  two  parameter 
Weibull  distribution.  It  should  be  noted  that  Lq  is  negative  as  A  and  S 
are  positive  numbers.  In  the  present  context,  this  parameter  does  not 
have  the  connotation  of  a  minimum  life.  It  is,  in  fact,  a  statement  of 
the  probability  of  a  first  cycle  failure.  In  particular,  if  N  •>  0 
Equation  41  reduces  to  the  static  distribution,  Equation  34 

Writing  Equation  42  in  the  form 


CSbN  = 
o 


f-mi 


with  C  ■  AS  ,  yields  the  basic  form  of  the  S-N  curve  which  is  different 
than  Equation  1.  It  should  be  noted,  however,  that  the  parameter  c  is 
simply  a  ratio  of  the  static  strength  shape  parameter  to  the  fatigue 
shape  parameter.  Experience  has  shown  this  number  to  be  typically  greater 
than  10  for  composites.  In  addition,  the  stress  levels,  S,  of  concern 
are  considerably  less  than  the  static  strength  location  parameter,  Sq. 
Thus, 


ft) 


and  Equation  44  becomes 


CS  N  =1 
o 


which  is  of  the  same  form  as  Equation  7.  Thus,  the  "wearout"  model 
approach  yields  an  S-N  curve,  for  stress  levels  of  interest,  which  is  of 
the  same  form  as  assumed  for  the  proposed  data  reduction  procedure. 

It  should  be  noted  that  the  "strength  degradation"  model  proposed 
by  Yang  and  Liu  (Reference  3)  is  not  based  on  a  flaw  growth  law.  Instead 
they  assume 

ASb=f(S,  wr)  <47) 


r 
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where  oi  is  the  frequency  and  r  the  stress  ratio.  They  are  forced  to 
assume  the  shape  of  the  S-N  curve  in  order  to  determine  f.  The  relation¬ 
ship  assumed  was  of  the  same  form  as  Equation  46.  Thus,  from  a  mathematical 
standpoint  their  model  is  the  same  as  the  crack  growth  or  "wearout"  model. 


Consider  Equation  39  in  the  form 


R(N)  =  exp 


where 


N  =A_1S  CS 
o  o 


b 


(48) 


(49) 


For  high  cycle  fatigue  the  second  term  in  the  exponential  is  negligible 
and  R(N)  can  be  approximated  by  the  two  parameter  Weibull  distribution 

R(N)  =  exp  (50) 


Thus,  the  data  reduction  procedure  proposed  measures  parameters  which 
are  consistent  with  "wearout"  models  or  "strength  degradation"  models. 

Such  models  must  be  modified  from  the  details  presented  here  if  more 
complex  loadings  other  than  T-T  are  to  be  considered.  Competing  failure 
modes  which  may  interact  are  a  problem  in  the  case  of  T-C,  for  example. 
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TABLE  3 

S-N  CURVE  AND  POOLED  WEIBULL  PARAMETERS 


Laminate 

b 

-- 

K/S 

o 

a 

£ 

X 

o 

l(T-T) 

66.34 

1.066 

0.310 

0.9490 

2(T-T) 

21.70 

1. 188 

1.08 

1.014 

2(T-C) 

9.705 

1.791 

1.45 

0.  9863 

3(T-T) 

65.84 

1.094 

3.82 

1.003 
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Figure  3.  Comparison  Between  Data  and  Pooled  Two  Parameter  Weibull 
Distribution  S/S  =  0.71 
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